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Abstract: We consider 'minimal' Z' models, whose phenomenology is controlled 
by only three parameters beyond the Standard Model ones: the Z' mass and two 
effective coupling constants. They encompass many popular models motivated by 
grand unification, as well as many arising in other theoretical contexts. This pa- 
rameterization takes also into account both mass and kinetic mixing effects, which 
we show to be sizable in some cases. After discussing the interplay between the 
bounds from electroweak precision tests and recent direct searches at the Tevatron, 
we extend our analysis to estimate the early LHC discovery potential. We consider a 
center-of-mass energy from 7 towards 10 TeV and an integrated luminosity from 50 
to several hundred pb _1 , taking all existing bounds into account. We find that the 
LHC will start exploring virgin land in parameter space for Mz> around 700 GeV, 
with lower masses still excluded by the Tevatron and higher masses still excluded by 
electroweak precision tests. Increasing the energy up to 10 TeV, the LHC will start 
probing a wider range of Z' masses and couplings, although several hundred pb" 1 
will be needed to explore the regions of couplings favored by grand unification and 
to overcome the Tevatron bounds in the mass region around 250 GeV. 
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1. Introduction and summary 

Among the extensions of the Standard Model (SM) at the TeV scale, those with an 
additional U(l) factor in the gauge group, associated with a heavy neutral gauge 
boson Z', have often been considered in direct and indirect searches for new physics, 
and in the studies of possible early discoveries at the LHC (for recent reviews and 
references, see e.g. [1]). While not prescribed by compelling theoretical or phe- 
nomenological arguments, these extensions naturally arise from Grand Unified The- 
ories (GUTs) based on groups of rank larger than four and from higher-dimensional 
constructions such as string compactifications. Z' bosons also appear in little Higgs 
models, composite Higgs models, technicolor models and other more or less plausible 
scenarios for physics at the Fermi scale. 

Many varieties of Z' models have been considered over the years [1]. In the 
following, we will concentrate on a class of minimal models, previously discussed in 
[2], that stands out for its simplicity and for the small number of additional free 
parameters with respect to the SM ones. Nevertheless, this class is sufficiently rich 
and motivated to emerge as a natural benchmark for comparing direct and indirect 
signals in different experimental contexts, in particular for organizing experimental 
searches at present hadron colliders such as the Tevatron and the LHC. 

By minimal Z' models we mean the most economical U(l) extensions of the SM 
that do not spoil renormalizability. Making reference to the SM particle content, our 
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minimality requirements can be summarized as follows: no exotic vectors, apart from 
a single Z' associated with a U (1) factor in the gauge group, commuting with G$m — 
SU (3)c x SU (2)lxU (l)y ; no exotic fermions, apart from one right-handed neutrino, 
singlet under Gsm, for each of the three SM families, and family- independent but 
generically non- vanishing U(l) charges; no exotic scalars beyond the SM Higgs field, 
meaning that the extra U(l) factor in G is either broken explicitly (which is still 
compatible with renormalizability) , or that the possible non-SM scalars from the 
extended Higgs sector can be neglected in the first Z' 'discovery studies', because 
they are sufficiently heavy and/or sufficiently decoupled from the SM fields, including 
the SM Higgs. 

The introduction of right-handed neutrinos makes it possible to generate a re- 
alistic pattern of neutrino masses (Dirac and Majorana) via renormalizable interac- 
tions. Similarly, it is possible to cancel all gauge and gravitational anomalies without 
introducing exotic fermions (as is instead required in some E 6 models) and/or non- 
renormalizable anomaly-canceling terms (as in some string-inspired models), as long 
as the generator of the additional U(l) factor is a linear combination of the weak 
hypercharge Y and of baryon-minus-lepton number, B — L. Also, our minimal class 
of models interpolates continuously among several discrete examples already consid- 
ered in the literature: the 'pure B — V model, the model arising from SO (10) 
unification, left-right symmetric models, etc. Minimal Z 1 models can be extended to 
include supersymmetry, but the price to pay is the introduction of many free param- 
eters associated with the supersymmetric particle spectrum: for this reason we will 
restrict ourselves to the minimal case, commenting on the possible effects of light 
supersymmetric particles when relevant. 

As we shall see, the simplicity of the parameterization allows us to perform a 
relatively model-independent study, which does not focus on a fixed type or size of 
Z 1 couplings and automatically takes into account mixing effects, both in the mass 
terms and in the kinetic terms, such as those generated in the evolution under the 
renormalization group equations (RGE). We find that the latter effects are indeed 
quite important, even for relatively weakly coupled Z' . 

With the parameterization above, we perform an updated analysis of the present 
indirect bounds coming from electroweak precision tests (EWPT), including LEP2 
and other experiments at low energy, and of the recent limits from direct searches at 
the Tevatron, taking into account the effects of mixing when necessary. Contrary to 
the common lore, we find that models with no Z — Z 1 mixing, such as pure B — L 
models, are as constrained as other models; actually, they are even more constrained 
than those with a partial cancellation of some suitable effective charge, such as the 
X model. 

Our goal is to apply our parameterization to study the impact of the present 
bounds on the discovery reach of the LHC, especially in the early phase when energy 
and luminosity will be limited. We find that the first virgin land in parameter space 
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that will become accessible to the LHC will correspond to relatively weakly coupled 
and light Z 7 , with masses of 600-800 GeV, while Z 1 bosons of the kind favored by 
GUTs, which are forced to be heavier by the present bounds, will require some 
more energy and luminosity to become accessible. We will show that indeed the 
unexplored region of parameter space — masses and types of Z' couplings — that will 
become accessible at each energy and luminosity, especially during the first LHC 
runs, depends non-trivially on the present bounds. This makes the use of a model- 
independent parameterization, such as the one suggested in this work, a valuable 
tool to systematically organize the Z 1 searches. 

Our paper is organized as follows. In Section 2, we elaborate on some results of 
[2] and illustrate how, in our minimal class of models, the main Z 1 properties and 
its couplings to the SM states can be very simply described in terms of just three 
new parameters beyond the SM ones: the Z' mass, Mz', and two effective coupling 
constants, gy and gsi, associated to the Y and B — L currents, respectively 1 . Such 
parameterization automatically takes into account all mixing effects in the gauge 
boson sector, both in the kinetic terms and in the mass matrix. In contrast with 
GUT-inspired parameterizations, often used for interpreting experimental searches, 
our parameterization also allows for relatively weak Z 1 couplings to the SM fermions: 
such a situation is disfavored by conventional GUTs but could arise, for example, 
in anisotropic string compactifications. To give a flavor of the generic GUT con- 
straints, we solve the one-loop RGE for the above coupling constants, both in the 
non-supersymmetric and in the supersymmetric case, showing that the large hier- 
archy between the GUT scale and the weak scale generically induces sizable kinetic 
mixing effects from the running, even for moderate values of the coupling constants. 
We then identify both the region of parameter space generically favored by GUTs 
and some points in parameter space corresponding to 'true' GUT-models, compar- 
ing them with some conventional benchmark points often used in the literature. In 
Section 3, we review the present bounds on the Z' parameter space. We start by 
discussing the impact of EWPT, including those from LEP2 and atomic parity viola- 
tion (APV), applying the methods and the results of [3]. In addition, we discuss the 
impact of some later reanalysis of APV [4], to conclude that the electroweak fit of [3] 
still gives the most stringent constraints. Then we discuss the bounds coming from 
direct searches at the Tevatron, including the most recent data from CDF [5, 6] and 
DO [7]. We confirm that, as for the models considered in previous studies [8, 9, 10], 
EWPT are more stringent than Tevatron direct searches for relatively heavy and rel- 
atively strongly coupled Z', such as the ones predicted by GUTs. On the other hand, 
Tevatron searches are more stringent than EWPT for relatively light and relatively 
weakly coupled Z', disfavored by GUTs but potentially permitted by other models. 
As a representative example, we comment on the possible excess in the CDF e + e~ 

1 The reader should not confuse gy with the SM U{l)y gauge coupling constant, which we denote 
by g', with g denoting the SU(2)l gauge coupling constant. 
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sample [5] at an invariant mass near 240 GeV: the size of the effect would correspond 
to a rather weakly coupled Z 7 , disfavored by GUTs but still allowed by EWPT. We 
conclude in Section 4 with a study of the prospects for the first LHC run(s). The 
main question we address is what energy and luminosity are needed to explore virgin 
land in parameter space, and how much of such accessible new territory is compati- 
ble with conventional GUTs. We show that, with the foreseen schedule for the first 
year [12] of the LHC (first 50-100 pb" 1 at ^ = 7 TeV, then up to 200^300 pb" 1 at 
< 10 TeV), the first region in parameter space to be explored will correspond to 
moderately light and weakly coupled Z', weighing around 600-800 GeV, and with a 
small window of allowed couplings. To open up considerably the region of parameter 
space accessible for discoveries, in particular the one relevant to GUT models, at 
least 0(1) fb^ 1 of integrated luminosity should be collected. In summary, in the 
very first phase of the LHC the interplay among center-of-mass energy, integrated 
luminosity and previous direct and indirect bounds will be quite subtle, and it will be 
important to focus the analysis onto the most promising regions of parameter space, 
possibly combining different channels and experiments from the very beginning. 

2. Theory 

2.1 Parameterization 

As discussed in the Introduction, we will consider extensions of the SM where Gsm 
is extended by a single additional non-anomalous family-independent U(l) factor, in 
the presence of three full SM fermion families, including right-handed neutrinos. As 
for the Higgs sector of the theory, we will assume the existence of the SM Higgs field 
but avoid as much as possible any specific assumption on the symmetry breaking 
mechanism for the additional U(l). 

As previously discussed in [2] , it is not restrictive to parameterize masses, kinetic 
mixing [13] and interactions with fermions for our extended neutral electroweak sector 
by means of the following effective Lagrangian: 

C = -\h AB F* F B ^ + \ M\ B A A »A* + A*r A + ..., (2.1) 

where A, B = T^l, Y, B — L. It is also well known that, by appropriate field redefi- 
nitions, we can go to a field basis where kinetic terms are canonical and masses are 
diagonal: 

C = -\f; v F*i» + \ M?A» 4 + 4 Jf + . . . , (2.2) 

where % — 7, Z, Z' . In the above equation, M 7 , Mz and Mz> are the mass eigenvalues. 
The currents Jf {% = 7, Z, Z') are those coupled to the gauge boson mass eigenstates. 
For example, 

J» = eJ2Q(f)frf (2.3) 
/ 
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Table 1: The charges of left-handed fermions controlling the electroweak neutral currents. 



is the electromagnetic current, where / runs over the different chiral projections of 
the SM fermions, Q(f) = T^,L{f)+Y{f) is their electric charge, and the contributions 
from the scalar sector have been omitted. Similarly, we can write 

J£ = cos & J£ - sin & J| /0 , J% = sin & J| + cos & J|, , (2.4) 

where 

J z» =9zJ2 PW/) - ^n 2 6w Q(f)] f-ff , (gz = , (2-5) 

/ 

is the SM expression for the current coupled to the SM Z° (we recall that, in the 
presence of mixing, Z° does not coincide with the mass eigenstate Z), and 

Jz'o = E & Y (f) + 9bl (B - L)(f)\ f-ff 

f 

= X><M/)/V7- (2-6) 
/ 

Again, possible contributions to the currents from the scalar sector have been omit- 
ted. We collected in Tab. 1 the charges of the SM fermions needed for evaluating the 
currents of eqs. (2.5) and (2.6). For definiteness, we chose a purely left-handed basis 
for the fermion fields, so that, omitting family indices, f — u,d, u c , d c , v, e, u c , e c . In 
expressing the charges Qz>, we found it convenient to make reference to the ratios 

~ 9y ~ 9bl / _x 

9y = j 9bl = • (2.7) 

9z 9z 



The parameterization above automatically contains and extends specific models 
often considered in the literature, such as Zb-l, Z x , and Z 3 b models, whose couplings 
simply read, in our notation: 
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where gz> is usually fixed to a 'GUT-inspired' value gz> = a/5/3 g'. 

Since the SM Higgs doublet H has 2 vanishing B — L, and, as discussed in [2], it 
is not restrictive to take the Higgs fields that break B — L (if any) to have vanishing 
Y, we can express the Z — Z' mixing angle 9' in terms of gy and Mz>, 



ta ' l(, '=-^MF%- (2 ' 9) 



where 

Mfo = &f (2.10) 

is the SM expression for the Z° mass. The same remains true if we assume that 
there is an explicit (or Stiickelberg-like) diagonal mass term for the Z ', without 
introducing an additional complex Higgs field for breaking B — L. 

Notice that the mixing angle is completely determined by the mass and the 
couplings of the Z'. In particular, it is always non-vanishing whenever gy ^ (i.e. 
for models different from pure B — L), because in these cases gauge invariance of 
the Yukawa terms forces the SM Higgs to be charged under the extra U(l), thus 
producing a Z — Z' mixing. 

We can then study the Z' phenomenology in terms of three unknown parameters: 
the Z' mass Mz> and the two coupling constants (gy,gBL) or, equivalently, (gy^Bi)- 

We will not consider possible additional parameters of the enlarged Higgs sec- 
tor and the right-handed neutrino masses, because, as will be discussed later, these 
parameters will play a relatively minor role in the following. To be definite, we will 
assume that there are three mostly left-handed neutrinos lighter than 0(1) eV and 
three mostly right-handed neutrinos heavier than Mz'/2, as in the see-saw mecha- 
nism, and that the physical components of the Higgs fields whose VEVs break B — L 
(if any) have negligible mixing with the SM Higgs and masses larger than Mz>- 

2.2 Constraints from grand unification 

One of the possible motivations for considering Z' models are GUTs, with or without 
supersymmetry. Through appropriate boundary conditions at the unification scale 
Mu and RGE on the running gauge coupling constants, GUTs can constrain the 
range of some low-energy Z' parameters, such as the coupling constants gy and gsi- 
The most stringent constraints can be obtained within specific models, where the 
full particle spectrum is specified and threshold and higher-loop corrections can be 
computed. Here, instead, we would like to remain as model-independent as possible 
within the general class of minimal Z' models. To this end, we will identify a GUT- 
favored region in the (gy,gBL) plane 3 according to the following procedure. 



2 This property is shared by the MSSM Higgs doublets Hi and H 2 - 

3 All the results are insensitive to the transformation (<?y, <7bl) — * {—gy, —9bl), without lack of 
generality we will thus consider only the upper half plane gBL > 0. 
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Table 2: The bAB coefficients for the four representative cases defined in the text. 



First, we recall that the gauge coupling constants g' , gy, gBL are related to the 
2x2 submatrix Hab with A,B = {Y, B — L} appearing in eq. (2.1), via 



h-AB 



/I 9y 1 \ 

(2.11) 



g' 2 9bl 9 ' 2 



9y 1 1 9y 1 



V 9bl9' 2 9 2 bl 9 2 bl9' 2 ) 
The matrix Hab obeys simple one-loop RGEs, which can be solved analytically: 



b AB , (M^ : 



h A B(M u ) = hAB(M z )--^- 2 \og^j , (2.12) 
where 

bAB = ^J2QfQf + s^ Q ^ ' (2 ' 13) 

/ * 

and s are the complex scalars in the theory. We collect some representative values of 
the bAB coefficients in Table 2, always including the contribution of the three fermion 
families of Table 1, and adding: 

(i) only the SM Higgs field, H ~ (1, 2, +1/2, 0); 

(ii) the SM Higgs field H plus a complex SM-singlet scalar ~ (1, 1, 0, +2); 

(iii) spin-0 superpartners for all the quarks and leptons in Tab. 1, spin-1/2 su- 
perpartners for all the gauge bosons, and two Higgs chiral superfields H\ ~ 
(1, 2, -1/2, 0) and H 2 ~ (1, 2, +1/2, 0), as in the MSSM; 

(iv) all the fields of the previous case, plus two extra chiral superfields 0i ~ 
(1,1,0,-2) and 2 ~ (1,1,0, +2). 

The quantum numbers in brackets denote, in a self-explanatory notation, the repre- 
sentations of SU(3) C , SU(2) L , U(1) Y and U(1) BL . 

For choosing the boundary conditions at the GUT scale Mu, we normalize all 
U(l) charges as in 5*0(10), and we take Mu = 10 16 GeV as a reference value. In 
typical GUTs, Mu can vary within approximately two decades around such refer- 
ence value, but the difference in our estimate of the GUT-favored region is of the 
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-0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 

Figure 1: GUT-favored region and some representative models in the (<?y , Qbl) plane, see 
the text for details. 

order of other threshold effects that we reabsorbed in the wide ranges we assume 
below for other parameters. Then, we compute the boundary value g'(Mu) using the 
phenomenological input g'(M z ) = e(M z )/ cos 8 W (M Z ), with a~^(M z ) ^ 128 and 
sin 2 6 W (M Z ) ~ 0.23, and the SM one-loop RGE. We then allow the Z' coupling at 
the unification scale au = 9u/{^) = 9z'(-^u)/(^), to vary within the generous 
bounds 

<«{/<—. (2.14) 

100 u 20 K J 

Taking into account that the SM RGE would predict au ~ 1/45, our upper and 
lower bounds leave a margin of more than a factor of two to account for threshold 
corrections, new particles at the TeV scale and other model-dependent effects. Cor- 
respondingly, we determine the GUT-favored region of the (gy^BL) plane by making 
use of the one-loop RGE of eqs. (2.11)-(2.13): the result is presented as the colored 
band in Fig. 1. The same figure also shows some dots that represent either some pop- 
ular GUT-inspired benchmark models considered in experimental analyses (the three 
empty dots and the corresponding dashed lines) or specific SUSY-GUT models with 
an extra U(l) (the three pairs of full dots). In particular, and in counter-clockwise 
order: the three dashed lines correspond to the three different models of eq. (2.8), 
when g Z i is left free to vary; the three empty dots correspond to the GUT-inspired 
normalization g z < = a/5/3 g'(Mz)- Instead the SUSY-GUT models are derived prop- 
erly, using the RGEs: they assume that the GUT group, say 50(10), is broken at 
My into the SM gauge group times an additional U(l) factor, with charges fixed as 
in eq. (2.8) at the GUT scale. For each of the three models (which correspond, in 
counter-clockwise order, to those in eq. (2.8)) we draw two black points, correspond- 
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ing to the RGE evolutions of case (iii) (outer points) and case (iv) (inner points). In 
both cases we get My ~ 2 x 10 16 GeV and ay ~ 1/24, within the bounds of eq. (2.14). 
We do not consider non-supersymmetric GUTs among our examples, because they 
would require the introduction of rather ad hoc exotic fields at intermediate mass 
scales to match the measured gauge coupling constants at the weak scale. Notice 
that the introduction of the extra chiral superfields 4>i an d 02 makes the coupling 
constants gy and Qbl more 'infrared free' but does not change ay, thus the values 
of gy and gsi at the weak scale are smaller. 

An important point to notice is that, even if we start form a 'pure B — U or 'pure 
T 3 r' model at My, the mixing effects in the RGE generate, through the resummed 
large logarithms, sizable corrections to the effective weak-scale couplings, as can 
be seen from the displacement of the black dots Z 1 ^^}^ from the corresponding 
dashed lines in Fig. 1. With enough running, a specific Z' at the unification scale 
can turn into a completely different one at the weak scale! Such effects make clear 
the advantage of considering the full (gy,gBL) plane to parameterize Z' searches. 
Notice finally that the direction in the (gy,gBL,) plane corresponding to the x model 
is quite stable under RGE, because running effects are caused by the MSSM Higgs 
superfields only. The RGE seem to exhibit an infrared attractor towards the region 
of parameter space with gy ~ — Ijbl- this fact will become even more interesting 
after compiling the experimental bounds. 



3. Present bounds 

3.1 Electroweak precision tests 

The use of electroweak precision tests to put bounds on the Z' parameters has a 
long history. Pre-LEP bounds [14] became much more stringent [15] after LEP1 
data at the Z peak, which strongly constrain the mixing angle 0'. However, also 
higher-energy LEP2 data and APV play a very important and complementary role 
in constraining the Z' mass for given couplings, as emphasized for example in [8, 3, 
9, 10, 11]. In [3] it was shown that the bounds from EWPT can be conveniently 
rewritten into bounds for the nine EW pseudo-observables (S, T, U, V, W, X, Y, 
5e g , 5C q ). For the case of extra Z' models, these pseudo-observables have simple 
expressions in terms of the Z' couplings and mass. Once such bounds are rewritten 
in terms of our parameters (Mz', gy, <?bl), for each Z 1 mass we can extract the 
corresponding 95% CL exclusion regions in the (gy,gBL) plane, as shown in Fig. 2. 
The region allowed by EWPT is the red one enclosed by each contour. Dashed lines 
and empty (full) dots remind us of the GUT-favored region and of the GUT-inspired 
(-derived) benchmark models. We remind the reader that EWPT only constrain the 
ratio gz'/M z >: this explains why the size of the red regions in Fig. 2 grows linearly 
with the Z' mass. 
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Figure 2: The regions of the {qy^bl) plane allowed by EWPT, at 95% CL, for M z > = 200, 
500, 1000, 1500, 2000, 2500, 3000 GeV (from inner to outer). The GUT-favored region is 
between the dashed lines. 

It is worth noticing that a recent theoretical re-analysis [4] of the most precise 
measurements on APV [16] was not included in the fit of ref. [3]. As will be discussed 
later, the new bounds from APV are significantly stronger than before, but not strong 
enough yet to compete with the result of the global fit. 

It is usually thought that pure B — L models are less constrained by EWPT 
because of the absence of Z — Z' mixing. Notice however that the region with gy — 
is not particularly favored, actually the region of parameter space least constrained 
by EWPT is that with gy — —Qbl- This feature can be understood by looking at 
the last row of Tab. 1, which shows that the Z' is less coupled to matter fields, thus 
less constrained by LEP2 bounds (and by Tevatron bounds as well, as we will see in 
the next section), roughly when g Y ~ —gsL- 

Notice also the correlation between the orientation of the GUT-favored region 
of Fig. 1 (between the dashed contours in Fig. 2) and the EWPT-allowed regions of 
Fig. 2. For all values of the couplings in the GUT-favored region, and in particular 
for the SUSY-GUT models represented by the full dots, the lower bound on Mz> 
is above 1 TeV. It is typically above 1.5 4- 2 TeV for the GUT-inspired benchmark 
models often considered in the experimental literature. For comparison, we report 
in Tab. 3 the bounds on the Z' masses for the particular choices of the couplings 
corresponding to the GUT-inspired and the SUSY-GUT benchmark points 4 of Fig. 1. 

4 The alert reader will notice that the bounds in Tab. 3 are numerically stronger than those from 
Fig. 2: this is simply because the figure refers to a 2-parameter fit (gy , <?bl), whilst the table refers 
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Table 3: 95% CL bounds on the Z 1 masses from EWPT, corresponding to the specific 
models represented by the nine points in Fig. (1). The Z'^ models are those represented 
by empty points, while Z'^ (Z 1 ^) corresponds to the three external (internal) black 
points; see the text for details on the choice of the effective couplings. 



Finally, we comment on the Higgs mass dependence of the bounds above. Since 
the SM Higgs mass is unknown, we may worry about the stability of our fits to EWPT 
with respect to varying the Higgs mass. Notice, however, that the dependence of the 
EWPT on the Higgs mass is only logarithmic and, although Z' bosons may help 
weakening the EWPT bounds on the Higgs mass [17], in our minimal models the 
preferred value for the Higgs mass is still below the LEP bound. Varying the Higgs 
mass within the 95% CL limit from EWPT (which corresponds to rrih ~ 200 GeV) 
produces only a tiny shift in the regions plotted in Fig. 2. To be definite, we chose 
m h = 120 GeV as a representative value. 

3.2 Tevatron direct searches 

Other important bounds on Z' parameters come from direct searches at hadron 
colliders, presently dominated by the Tevatron experiments CDF and DO. In these 
experiments, Z' bosons of sufficiently low mass can be produced on-shell and decay 
in the process qq — > Z' — > £ + £~, (£ = e, //). These two are very clean channels 
to look for: a peak in the invariant mass distribution of e + e~ or pairs, with 

a width controlled by the experimental resolution when the intrinsic Z 1 width is 
sufficiently small. The irreducible background is dominated by SM Drell-Yan (DY) 
£ + £~ production, which is well understood and whose control is only limited by PDF 
uncertainties. Other irreducible backgrounds are small and the reducible ones can 
be easily eliminated by generous cuts. Different decay channels into SM final states 
have been also experimentally investigated [18], for example T + T~,jet jet,W + W~ , 
but for the minimal models considered here they are not competitive for exclusion 
or discovery: at best, they could play a role in the determination of the Z' couplings 
after a future discovery. 

The procedure used to extract bounds from such processes has been extensively 
discussed in the literature (see for example [9, 19]). For a weakly coupled Z', it 
consists in calculating the Z 1 production cross-section multiplied by the branching 
ratio into two charged leptons a(pp — > Z'X) x BR(Z' — > £ + £~) (in our case a function 

to specific models, thus to a 1-parameter (M' z ) fit. 
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of the three Z' parameters M z >, gy, 9bl), and comparing it with the limits established 
by the experiments. 

On the theory side, we performed the calculation at NLO in QCD (and LO for 
the EW part), using the NLO MSTW08 PDF sets [20]. In the calculation of the total 
width Y Z ' we included the following channels: Z' — > //, W + W~, and Zh, where h is 
the SM Higgs boson and / are the SM fermions of Tab. 1, with the exception of the 
right-handed neutrinos, which we took to be heavier than M z >/2. The presence of 
the two last decay channels is due to Z — Z' mixing and is usually neglected, however 
for large Z 1 masses there is an enhancement that cancels the suppression due to the 
mixing [21, 22]. The ratio Y Z i jMz> is pretty constant over the whole range of masses 
of interest, and is around 2% for GUT-favored Z' couplings, and of course smaller 
for more weakly coupled Z' . 

Notice that in the presence of extra matter fields charged under the extra U(l) 
(light right-handed neutrinos, more light Higgses, supersymmetric partners, etc.), 
Y z > would be larger, with a consequent suppression in the branching ratio to charged 
leptons. In this case the bounds from hadron colliders (and their ability for discovery) 
would be weaker (unlike those from EWPT, which are quite insensitive to these 
extensions of the model). 

For the Tevatron experimental limit we used the most recent available results 
from CDF (on Z' -> e+e" [5] and Z' -> [6]) and DO (on Z' -> e+e" [7]). They 

directly provide the 95% CL bounds on the product a(pp — > Z' X) x BR(Z' — > 
based on 2.5, 2.3, 3.6 fb _1 of data with 27^-38%, 13-j-40%, 17-^22% total acceptances 
respectively, with the acceptances growing from smaller to larger values of Mz>- 
Notice that, although DO data refer to a higher integrated luminosity, the acceptance 
is smaller, making its bounds a little weaker than those from CDF. 

Since not enough information is available to us to properly combine the three 
sets of data, for each Z 1 mass we took the strongest bound among the three sets. 
A combined analysis would be highly welcome as it would probably give stronger 
bounds. 

When compared with the computed cross-section, the experimental limits pro- 
duce, for each value of M Z ', a 95% CL exclusion region in the (<?y, 9~bl) plane, in 
analogy with the EWPT case. The results are summarized in Fig. 3, which shows 
the 95% CL allowed regions for different M' z . Notice that the region favored by 
GUT models starts becoming accessible for M' z > 700 GeV. For masses larger than 
~ 1.2 TeV the allowed region fills the whole plot, and the available data are no longer 
able to give useful constraints. Indeed, unlike the EWPT case, the allowed regions in 
Fig. 3 grow faster than linearly with M Z i, because of the stronger suppression from 
the x-dependence of the PDF at higher energies. 

3.3 Comparison among different bounds 

We may wonder whether it makes sense to compare direct and indirect experimental 
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Figure 3: The regions on the (flY, 9bl) plane allowed by Tevatron direct searches at 95% 
CL for M z > = 200, 400, 600, 800, 1000 GeV (from inner to outer). The GUT-favored 
region is between the dashed lines. 

bounds from different experiments, or better whether it is possible to build models 
that can evade indirect bounds but still be accessible to direct searches. Notice how- 
ever that indirect LEP searches, low-energy APV experiments, direct and indirect 
searches at hadron colliders are all basically controlled by tree-level Feynman dia- 
grams built from two basic types of elementary vertices, coupling the Z' to charged 
leptons and quarks, respectively: 



e e 




Of course, LEP and APV probe off-shell Z' exchange, whereas the Tevatron and the 
LHC are sensitive to on-shell Z' production and decay. But the parameters involved 
in the relevant Feynman diagrams are the same, and it is not easy at all to invent 
new physics capable of evading indirect bounds from EWPT but still producing a 
signal in the direct searches. For this reason, the bounds from EWPT should not 
be neglected when analyzing the discovery reach of direct searches. On the other 
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hand, if the branching ratio to leptons is suppressed by the presence of extra charged 
matter, then it may happen that indirect searches become even more powerful than 
the direct ones. 

The first thing we notice by comparing the plots in Figs. 2 and 3 is that both 
EWPT and Tevatron bounds probe similar regions of the (gy,gBL) plane. In partic- 
ular, they are both less stringent in the region 5 with g Y ~ ~9bl, and more stringent 
in the orthogonal direction (although the actual shape of the exclusion region from 
the Tevatron is slightly different from that of EWPT). As explained in the previ- 
ous section, this correlation in the couplings is due to the fact that in the region 
g Y ~ —qbl there is a partial cancellation in the Z'-charges of the SM fermions (see 
Tab. 1). Interestingly, also the region preferred by GUT models has the same shape 
as those allowed by the experimental bounds: this makes the bounds on M z > at fixed 
coupling only weakly dependent on the specific nature of the Z'. 

Another important difference between the bounds from the Tevatron and from 
EWPT is that, as already noticed in some previous analyses [8, 3, 9, 10], the former 
give stronger constraints at smaller M z > but weaker at larger M z >- Indeed, the 
possibility to produce on-shell Z'-bosons rewards hadron colliders in the mass region 
that is easily accessible to them. On the other hand, their power rapidly falls off 
at higher masses because of the PDF dumping at high x. Already for M z > around 
800 GeV, the EWPT bounds start becoming more powerful than the Tevatron limits, 
independently of the Z' couplings. In particular, all the models with GUT-favored 
couplings are bound much more strongly by EWPT than by direct searches, at least 
with the current data. 

The competition between the Tevatron and EWPT is more manifest if we plot 
the bounds as functions of M z * for a representative model. Fig. 4 illustrates the range 
of Z' masses where the Tevatron gives stronger or weaker constraints compared to 
EWPT, for a particular direction in the (^y , ~Qbl) plane corresponding to the so- 
called Z x models (the central dashed line in Fig. 1). Notice that for GUT-favored 
couplings EWPT give much stronger bounds than the Tevatron. Similar figures were 
previously shown in refs. [8, 10]. 

In Fig. 5, as an illustration, we combine the bounds from EWPT, Tevatron 
direct searches and the new APV analysis of [4], for the two representative values 
M Z i = 400 and 800 GeV. In the first case the bounds from Tevatron are the strongest 
for all models, while the APV bounds are the weakest. In the second case we see that 
EWPT start becoming stronger than those from the Tevatron in almost all parameter 
space: for example, Tevatron would still allow an 800 GeV Z' 3R with GUT-favored 
coupling, while EWPT would basically rule out all GUT-like Z' with this mass. 
APV is always weaker than EWPT but starts becoming stronger than the Tevatron 
in the region where the latter starts performing worse, i.e. the GUT-preferred region. 

5 Rcmarkably this region seems also to be an attractor of the solutions to the RGE, as commented 
in the previous section. 
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Figure 4: The region of the (Mz',9z') plane excluded at 95% CL by EWPT (red) and 
Tevatron (blue), for the Z x model defined in eq. (2.8). The horizontal strip (yellow) recalls 
the GUT-preferred region and the dotted line the value gz> = \/5/3 g'(Mz) corresponding 
to the Z'^ model. 




Figure 5: The region of the {gY,9Bh) plane allowed by EWPT (red), Tevatron (blue), 
and APV (green), for M z , = 400 (left) and 800 GeV (right). 



Notice also that APV experiments are not sensitive to certain types of Z' ', such as 
pure B — L that has purely vectorial couplings to fermions. However, as discussed 
in section 2.2, and evident in Fig. 5, RGE effects can easily move the Z' couplings 
away from the safe region, making such Z' models, which would be otherwise safe 
with respect to APV bounds, also subject to constraints. 

The CDF e + e~ excess 

Recently the CDF collaboration [5] observed a small excess of e + e~ events around 
240 GeV, amounting to a 2.5 a fluctuation from the SM background. However, CDF 
did not see any anomaly in the dimuon spectrum [6]. DO data [7] at the moment do 
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not seem to confirm nor exclude such excess, also because of the smaller acceptance 
mentioned before. Although our minimal models would not be able to explain such 
effect, because of the mismatch between electron and muon spectra, we may still 
ask whether existing bounds from EWPT may rule out an explanation of the excess 
in terms of some more general Z' bosons that couple non-universally to leptons, 
such as those discussed in [23]. Notice from Fig. 4 that at 240 GeV Tevatron is 
indeed more sensitive than EWPT to neutral resonances, provided that they have 
a small coupling, of order gz> ~ 0.04. It turns out that such small coupling would 
be enough to explain the Tevatron excess without contradicting EWPT 6 . We have 
not checked whether such scenario is compatible with other flavor non-universal low- 
energy constraints. More data are anyway required to confirm the presence of a true 
excess and assess its possible non-standard origin. Notice instead that an analogous 
signal at energies higher than ~700 GeV would not be compatible with EWPT, at 
least within the class of models considered here. 

4. Early LHC prospects 

The questions we want to address in this section are the following. At what combined 
values of center-of-mass energy and integrated luminosity may we expect the LHC 
to start having a chance of discovering a Z' (at least of the kind discussed in this 
work), taking into account all the experimental bounds discussed in the previous 
sections? What region of parameter space that has not been already ruled out 
could be accessible for different luminosities and energies in the first LHC runs? 
Considering the fact that Z' signals are among the cleanest and easiest ones in the 
search of new physics, our analysis may also be used as a benchmark point when 
discussing the integrated luminosities that are worth collecting at each energy to 
actually probe new physics. 

At the moment, the program for the first year of LHC running [12] consists in 
a very first run at low energy (y/s = 7TeV)) and low luminosity (< 100 pb _1 ), 
followed by an upgrade in energy (y/s < lOTeV), with a collected luminosity up to 
200^300 pb- 1 . 

At such low 7 energies and luminosities, the constraints from Tevatron direct 
searches and EWPT play a crucial role in identifying the unexcluded region of pa- 
rameter space that can be probed and the time scale required to have access to 
it. 

Being the LHC a hadron collider, the region of parameter space accessible to 
it will be similar in shape to the corresponding one at the Tevatron. For relatively 
light Z 1 (Mz> < 800 GeV), since the strongest constraints come from Tevatron direct 

6 Notice that the same statement can be true also for Z' models very different from those studied 
in this work [24]. 

7 Of course, with respect to the LHC design parameters. 
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Figure 6: The geometrical acceptance for signal (solid lines) and SM-DY background 
(dashed lines), as a function of a parameter that scans over the minimal models, and for 
two representative values of M^+£-: 200 GeV (red, lower) and 1 TeV (blue, upper). The 
different lines refer to the cut \rj\ < 2.5 and pxe > 20 GeV (thin) or p?i > 80 GeV (thick). 
The colored bands show how much the acceptance varies by changing the rapidity-cut from 
\r}\ < 2.1 to \r/\ < 3.0. 



searches, we expect the LHC to turn into a discovery machine as soon as it becomes 
sensitive to regions of parameters not yet excluded by the Tevatron. However, while 
the higher energy is clearly a big advantage for intermediate Z' masses of several 
hundreds GeV, for lighter masses the low luminosity may be a crucial limiting factor 
in the early LHC phase. On the other hand, for heavier Z' masses, such as those 
relevant for GUT models, generically EWPT outperform the Tevatron, and the LHC 
must wait for higher energies and luminosities to become sensitive. 

To turn these considerations into more quantitative statements, we perform a 
basic analysis along the lines of the one described before for extracting the Tevatron 
bounds. In the present case we consider the range i/s = 7 -j- 10 TeV for the pp 
center-of-mass energy, luminosities in the range 50 pb _1 4- 1 fb _1 , and calculate the 
product a(pp -> Z'X) x BR(Z' -> i+f) for M z > = 200 3000 GeV, at the same 
order in perturbation theory as in the Tevatron case. At the same level of precision, 
we also compute the SM Drell-Yan (DY) differential cross-section, which constitutes 
the main source of background. 

To gain some approximate understanding of the acceptances for signal and back- 
ground at different values of the invariant mass Mi+£- of the £ + £~ pair, and of the 
possible model-dependence of the former, we performed a simple study whose re- 
sults are illustrated in Fig. 6. For a/s = 10 TeV, we plot the purely geometrical 
acceptance for signal (solid lines) and SM-DY background (dashed lines), impos- 
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ing the cut \r)\ < 2.5. The colored bands show how much the acceptance varies if 
the cut is varied from \r)\ < 2.1 to \r]\ < 3.0. The quantity on the horizontal axis, 
[ArcTan [gy j Qbl)\I ^ scans over the different minimal models. The upper blue lines 
are for Mg+g- = 1 TeV, the lower red lines are for Mg+g- = 200 GeV: as expected, 
the acceptance for both signal and background is similar, grows with the invariant 
mass of the lepton-antilepton pair, and is close to 90% for Mg+g- of order 1 TeV 
or larger. We also looked at how the acceptance depends on the cut on the lepton 
transverse momenta pre'- the thin lines correspond to pTg > 20 GeV, the thick ones 
to pre > 80 GeV. We can see from Fig. 6 that this cut is essentially included in the 
cut on t] for high mass values, whereas it has a small but non-negligible effect on 
the acceptance for low mass values. Finally, the model-dependence also decreases 
when moving from lower to higher masses, and is never larger than 10% even for 
M z > ~ 200 GeV. In view of these results, and for the purposes of the present ex- 
ploratory study, we then assumed an acceptance depending only on the invariant 
mass, as done for example in ref. [27]. Our computed values of the acceptance 
are compatible with those of ref. [27], thus we adopt their Fig. 2 for the rest of 
our LHC study. More refined studies, however, should take into account also the 
model-dependence of the acceptance, which may not be negligible for Z' searches at 
relatively small masses. 

To estimate the 5a discovery reach of the early phase of the LHC [25, 26, 27], 
we compared the events due to a generic Z' signal to the events from the SM-DY 
background in a 3% interval around the relevant values of the dilepton invariant 
mass 8 . We then required the signal events to be at least a 5a fluctuation over the 
expected background, and in any case more than 3. This rough statistical analysis is 
enough to get an approximate answer to the questions we want to address. We leave 
a more careful analysis to the experimental collaborations ATLAS and CMS, which 
have control on all the information needed to perform it in an accurate and reliable 
way. A more refined analysis would also be needed for a possible Z' diagnostics after 
discovery, as studied for example in [28]. 

We present sample plots in Figs. 7 and 8. Besides the regions allowed by EWPT 
and Tevatron data (red and blue) we plotted, for each representative value of the Z' 
mass, and of the LHC energy and integrated luminosity, the region not accessible to 
the LHC (in yellow) for a 5a discovery as defined above. We see that \fs = 7 TeV and 
50 pb _1 are not enough to discover any Z' in the whole parameter space considered; 
in particular at low masses the low luminosity makes the LHC underperform with 
respect to the Tevatron, while at masses where the LHC starts having a kinematical 
advantage over the Tevatron, both cannot compete anymore with EWPT. With 
100 pb _1 of data at the same energy, a first non-excluded region of parameters 
becomes accessible to discovery, though it is very narrow (M Z / ~ 700 ± 100 GeV and 

8 This should be compatible with the expected energy resolution, even in this early phase, and 
with the fact that, for GUT-favored values of the coupling constants, Tz'/Mz' ~ 2%. 
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Figure 7: The LHC 5a discovery potential in the (gy,9BL) plane for y/s = 7 TeV. The red 
and blue regions are those allowed by EWPT and Tevatron bounds respectively; the yellow 
region is the one not within 5a discovery reach at the LHC. Thus the region accessible by 
the LHC is the one formed by points that are both in the red and blue regions but not 
in the yellow one. Plots in the first row refer to 50 pb _1 of data and Mz 1 = 200, 500, 
700 GeV respectively; plots in the second row are for 100 pb" 1 of data and M z > = 600, 
700, 800 GeV respectively. 




Figure 8: The LHC 5a discovery potential in the (<?y, 9bl) plane for y/s = 10 TeV, 
200 pb- 1 of data and M z , = 400, 500, 900, 1300, 1400, 1600 GeV. The meaning of the 
colored regions is as in Fig. 7. In the last three plots the Tevatron bounds are not shown 
because they are too weak to give useful constraints. 

g BL ~ 0.15 -=- 0.20, g Y ~ -0.2 0). 

Things start improving as the LHC steps up in energy and luminosity. The 
situation with \/s = 10 TeV and 200 pb -1 of integrated luminosity is represented 
in Fig. 8. The region of Z' masses below 400 GeV will not be accessible yet, this 
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M r (GeV) Mr (GeV) 

Figure 9: First row. The region of the (Mz',gz') plane amenable to a '5<r' discovery at 
the LHC, for the Z x model, y/s = lOTeV and some representative values of the integrated 
luminosity; from left to right: 50, 100, 200, 400 and 1000 pb -1 . The red and blue region 
and the yellow band are the same as in Fig. 4. The second box is a zoom on the low- mass, 
low-coupling region. Second row. 95% CL exclusion contours from the LHC after 50 and 
100 pb" 1 at y/s = 7 TeV (blue curves) and after 50, 100 and 200 pb" 1 at y/s~ = 10 TeV 
(green curves). 

because the higher luminosity collected at the Tevatron is more important in such 
energy region. The first accessible zone in parameter space starts showing up for 
Mz> ~ 400 -T- 1100 GeV, for models with couplings smaller than those preferred by 
GUTs, and for M z > ~ 1200 1500 GeV for GUT-like couplings. For heavier Z' no 
region is left to the LHC that is not already ruled out by EWPT. As evident from 
the plots, for each of the accessible Z' masses, only a small portion of the (gy^BL) 
plane will be tested. Our plots refer to data collected by a single experiment and 
for a single dilepton channel, combining the data might help increasing the effective 
luminosity collected and thus the discovery potential. 

Since contour plots may require some patience to be interpreted, we make our 
results more manifest by plotting, in Fig. 9, the 5a LHC discovery potential in the 
(Mz'iQz') plane for the representative \ model. As in Fig. 4, the red and blue regions 
are those presently excluded by EWPT and Tevatron direct searches, respectively, 
and the yellow band denotes the GUT-favored region. The new curves enclose the 
region where a 5a discovery at the LHC is in principle possible, for y/s = 10 TeV 
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and some representative values of the integrated luminosity: 50, 100, 200, 400 and 
1000 ptT 1 , from left to right. Notice that, in the case under consideration, the first 
mass region to be touched is between 600 and 800 GeV, with the region enlarging 
towards higher mass values with increasing luminosity. We start entering the GUT- 
favored region of parameters only for 100 pb _1 , and 1 fb _1 is enough to reach mass 
values as high as 2 TeV, with a full coverage of the GUT-favored region of couplings. 
Notice also that the access to lower mass values is also gradual, and that 1 ftr 1 is 
required to do better than the Tevatron at Mz> ~ 200 GeV. We have performed a 
similar analysis for y^s = 7 TeV: in such a case, 400 pb -1 give approximately the 
same sensitivity as 200 pb -1 at 10 TeV for Mz> < 700 GeV, whilst the sensitivity 
at higher mass values rapidly becomes worse and worse, as expected: there are no 
doubts that it is worth raising the LHC energy as soon as it can be safely done. 

Of course, if no discrepancy from the SM is found in the dilepton spectra, LHC 
will be able to improve the 95% CL bounds on minimal Z' models already after the 
first run(s). Indeed, as shown in Fig. 9, 100 pb~ x at y/s = 7 TeV are already enough 
to top both Tevatron and EWPT bounds for all Z' masses up to 1.3 TeV, while, 
after the first year of run, LHC might be able to rule out most of the GUT-preferred 
region below ~ 2 TeV. 

In summary, our study shows how strong the role played by the existing experi- 
mental bounds can be in limiting the access to new physics in the early LHC phase. 
Even for the 'easy' Z 1 models usually considered, the energy and luminosity required 
to overcome existing bounds can delay the possibility of discovering new physics by 
a non- negligible amount. The importance of reaching higher energies and luminosi- 
ties is clear, as it is the importance of combining data from different detectors and 
channels already in this early phase. 

Our analysis also shows that different regions in parameter space will become 
available for discovery at different times, depending on the energies and luminosity 
reached. Hence it would be sensible to switch to general parameterizations such as 
the one described in the present work, which are not as restrictive as those commonly 
used at the moment. The latter may focus the attention on the 'wrong' regions of 
parameter space that are already ruled out by current data, instead of those with 
the greatest potential to host accessible new physics. 
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